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NATIONAL ADVTSOET COMMITTEE FOR AERONAUTICS 
KACA TN No. 1680 


DIVERGENCE OF SWEPT WINGS 
By Franklin W. Diederich and Bernard Budiansky 


SUMMAET 


An analysis of the divergence of swept untapered and tapered wings 
with stiffnesses varying as the fourth power of the chord has "been 
performed and checked experimentally. The results are presented in a 
set of charts and approximate formulas suitable for quick estimates of 
the divergence dynamic pressure and hence the divergence speed. 

These results indicate that the divergence speed drops rapidly 
as sweepforward increases to about 40° but that wings with moderate or 
large sweepback cannot diverge. The location of the elastic axis is 
found to affect the divergence speed most at low angles of sweep, where 
movement of the elastic axis forward (or the aerodynamic center aft) 
haises the divergence speed. The effect of wing taper is to increase 
the divergence speed of essentially unswept wings and to decrease the 
divergence speed of wings with moderate and large angles of sweep in 
the case of the prescribed stiffness variation. Evidence is presented 
to indicate that these effects may not be ^observed for actual stiffness 
variations, in which cases a more refined analysis must be resorted to. 


INTRODUCTION 


The- emphasis on the use of sweptback or sweptforward wings for 
high-speed flight has created widespread interest in the aeroelastic 
behavior of swept wings. The present paper is concerned with the 
theoretical determination of one of the most fundamental aeroelastic 
parameters, the wing divergence speed. 

The divergence of wings or tail surfaces is an instability 
phenomenon which results from the interaction of aerodynamic and 
structural forces. If a wing or tail is given a deflection of 
arbitrary magnitude, the aerodynamic forces often act in such a way 
as to increase the given deflection, whereas the structural forces 
always tend to decrease the deflection. Since the aerodynamic forces 
increase with the flying speed, whereas the structural forces are 
independent of it, a speed will often exist at which the two sets of 
forces are exactly in balance, so that they tend to maintain the given 
deflection. This speed is known as the divergence speed, since any 
further increase in speed causes the aerodynamic forces to predominate 
over the restraining structural forces and tends to increase any 
deformation until structural failure occurs. 
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The theory of divergence of unswept wings has reached a considerable 
degree of refinement. Since H. Beissner's original analysis of the 
problem, (reference l), a great deal of material has been published on 
this subject both in the United States and in Great Britain (for instance, 
references 2, 3j and 4); the latest methods treat arbitrary stiffness 
variation (reference 2) and account for aerodynamic induction 
(references 3 and 4} . 

The analysis of the divergence of swept wings is complicated by 
the fact that, unlike the case of unswept wings, the air forces depend 
on the bending deformations as well as on the twisting deformations. 

A first approximation to the solution of the problem is presented in 
reference 5 "by means of a "semirigid" approach that does not take into 
account spanwise variation of the wing distortion. 

The present paper gives the results of a more exact analysis that 
considers the effect of elastic bending and twisting along the wing span 
Theoretical derivations are given for the divergence speeds of swept 
wings with constant chord and constant flexural and torsional stiffnesses 
along the span, as well as for swept wings with linearly varying chords 
and with flexural and torsional stiffnesses varying quartically with 
the chord. The results of the theoretical analysis are presented in 
curves of nondimens ional parameters from which the divergence speed can 
be estimated for a given design. 

In order to verify the theory and establish the effects of the 
assumptions involved, a limited number of tests were made in the 
Langley 4.5— foot flutter research tunnel' on models of constant chord 
and stiffness at low Mach numbers. The results of these tests are 
presented and are compared with the theory. 


SYMBOLS 


A 


D 


aspect ratio 


Spanf 


k Total area 


dimensionless parameter 


cos%m e e-[_ c r 2 L 2 


v, (GJ) r 

value of parameter a at divergence 

cos‘v\m e c r L3tan 


dimensionless parameter 


(El). 


value of parameter d at divergence 

chord, measured perpendicular to elastic' axis, feet 
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'•'r 

El 


(El) ri 

®1 

GJ 

(GJ) r 

k 

Ki, K 2 
2 
L 

M ■ 

Mcr 

“o 

“e 

a 


V D 


chord at effective root (fig. l), feet 

effective chord at tip (fig. l), feet 

effective "bending stiffness in planes perpendicular to 
elastic axis* pound— feet 2 

bending stiffness at effective root, pound— feet 2 

distance from elastic axis to aerodynamic center 
(positive forward), fraction of chord 

effective torsional stiffness in planes perpendicular to 
elastic axis, pound— feet 2 

torsional stiffness at effective root, pound-feet^ 

chord ratio, c/o^. 

constants 

running air load along elastic axis, pounds per foot 
length of one wing along elastic axis (see fig. 1), feet 
free— stream Mach number 

critical Mach number of section perpendicular to 
elastic axis 

section lifi>-curve slope, per radian 
effective section lift-curve slope, per radian 
dynamic pressure, pounds per square foot 
dynamic pressure at divergence, pounds per square foot 

free— stream velocity, feet per second 

free— stream velocity at divergence, feet per second 

component of free— stream, velocity normal to elastic axis, 
feet per second 

distance along elastic axis (see fig. l), feet 
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T) 

r 

A 

x 

p 

<p 


effective angle of attack of a section perpendicular to 
elastic axis, radians 

semispan position along elastic axis (V/k) 

local dihedral angle* radians * or slope of -wing deflection 
cinpre at elastic axis (see fig. l) 


angle of sweep measured to elastic axis* positive for 
sweepback* degrees 


taper ratio 



density of air stream* slugs per cubic foot 


angle of twist in a plane perpendicular to elastic axis* 
radians 


THEORETICAL RESULTS 

f 

Untapered Swept Vings 


A theoretical analysis of the divergence of a swept wing of uniform 
chord (fig. 1(a)) and stiffness is contained in the appendix . The 
analysis involves the following limitations and assumptions: 

(a) Aerodynamic induction is taken into account only insofar as 
an over-all correction is applied to strip theory 

(b) Aerodynamic as well as elastic forces are based on the 
assumption of small deflections 

(c) The wing is clamped at the root perpendicular to a straight 
elastic axis (see fig. l(a) ) * and all deformations are given by the 
elementary theories of bending and of torsion about the elastic axis . 


Within the limitations of these assumptions an exact theoretical 
solution for the divergence speed is obtained. The solution consists 
of a relationship between two nondimensional parameters* 




o p p 

q. D cos Amge^Cj, 

Wh 


( 1 ) 
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and 



< 3 p cos^Am e Cpl.3 tan A 

(El) r 


( 2 ) 


■which is presented in figure 2. 


It is seen that the ratio of these parameters 


*D 

a D 


(GJ) r " 

L “ 

_(ED r _ 

e l c r 


tan A 


(3) 


depends only on the known geometrical and physical parameters of the 
problem. Thus, from the theoretical plot of ap against dp/ap shown 

in figure 3, the divergence speed can he obtained for any particular 
uniform wing. In order to cover the entire range of values of the 
independent parameter, it is convenient to plot dp against ap/dp for 

large positive and negative values of dp/ap. 


Tapered Swept Wings 

A theoretical solution for the divergence speed of linearly tapered 
swept wings (see fig. 1(b)) has also been effected in the appendix. 

The assumptions listed in the foregoing section apply to this case as 
well. The bending and torsional stiffnesses are assumed to vary as 
the fourth power of the chord; this variation is realized for a wing 
having geometrically similar cross sections, such as a solid wing or one 
all structural dimensions of which are proportional to the local chord. 

Theoretical curves similar to those of figure 3, relating the 
nondimens ional parameters ap and dp, have been computed for swept 

wings having taper ratios X of 0 . 2 , 0 * 5 , and 1 . 5 * The various branches 
of these curves, together with those for the .uniform wing (taper 
ratio 1.0) are, for convenience, given separately in figure 4(a), 4(b), 
and 4(c) . 


Effective Lift-Curve Slope 

In the calculations of the lift acting on- the individual wing sections, 
it has been assumed that the aerodynamic interaction of the sections nay be 
neglected if an over— all span correction is made to the section lift-curve 
3lope. Thus, the value of m e to be used' is the two-dimensional value Hq 
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for sections normal to the elastic axis multiplied by a suitable correction 
factor. In reference 2, Shornick, in computing the divergence speed of 
unswept wings by strip theory, applies the aspect— ratio correction 


“e 



(4) 


For the case of an unswept uniform wing with ~ = 3.14, this assumption 

yields a divergence speed that differs by less than 1 percent from the 
result calculated by Hildebrand and Beissner (reference 3) by lifting— line 
theory. (The parameter 0 calculated by lifting— line theory in reference 3 


may be compared directly with \J~&q 



of the present paper.) 


Equation (4) may be extended to swept wings by means of the reasoning 
used in reference 6. If the air forces are considered to act in planes 
perpendicular to the elastic axis instead of parallel to the plane of 
symmetry, the following relation is obtained 


®e - “o ~ (4a) 

A + 4 cos A 

This approximation for finite— span effect may be used for subsonic and 
subcritical Mach numbers. For supersonic and supercritical Mach numbers, 
it is inapplicable; no span correction is available at present for these 
speed ranges. Neglect of the finite— span effects, however, will always 
tend to give conservative results . 


Location of the Effective Boot 

If the results of the preceding analysis are applied to an actual 
wing, it is necessary to assume that the wing is clamped along a line 
perpendicular to the elastic axis. From the data and the analysis 
presented in reference 7 it appears that the amount of wing twist can 
be estimated closely by assuming the wing to be clamped at a line 
‘-trough the intersection of the elastic axis and the side of the fuselage, 

1? the carry— through structure is fairly rigid. 

The bending deflections are estimated in reference 7 by considering 
— ie wing supported flexibly at a line through the innermost point of the 
'- ra -ling edge (in the case of a swepthack wing) or of the leading edge 
( rn the case of a sweptforward wing) . This effect may be taken into 
account in the preceding analysis by modifying the root boundary conditions. 
If the carry— through structure is fairly rigid the resulting divergence 
speed is the same, for all practical purposes, as the value obtained by 
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considering the wing mounted rigidly at a line through the intersection of 
the elastic axis and the side of the fuselage. This line is consequently 
considered the effective root of an actual wing, as shown in figure l(c). 


APPLICATION OF THE THEORETICAL RESULTS 
Selection of the Aerodynamic Parameters 

f 

Both the section lift— curve slope hiq • and the aerodynamic center 
(and consequently the distance from the aerodynamic center to the elastic 
axis, ep) vary with Mach number and, in the transonic range at leasts 
with the angle of attack. The parameters should he chosen at the angle 
of attack near the design value which yields the most conservative results. 
The choice of the Mach number depends on the purpose for which the 
divergence speed or dy nami c pressure is calculated. 

If the divergence speed is calculated for its own sake, the parameters 
should be chosen at the Mach number corresponding to the divergence speed; 
resort to a trial and error procedure may therefore be necessary. If, on 
the other hand, it is desired to calculate the divergence speed or dynamic 
pressure as a reference value for one of the other aeroelastic phenomena, 
the aerodynamic parameters should be chosen at the Mach number of interest 
for the particular phenomenon. The dynamic pressure calculated in this 
manner will not usually be the true value but appears to be the divergence 
value only at the given Mach number and will vary with Mach number. If 
the variation of this reference divergence dynamic pressure qp is 
plotted against Mach number and the actual dynamic pressure q is plotted 
on the same graph, the intersections of the two lines will determine the 
true values of the divergence Mach number and dynamic pressure; this 
procedure therefore constitutes another way of calculating the divergence 
speed.- Such a plot is shown in figure 5 for a straight and a sweptforvard 
wing designed for highr-speed flight. 

The values of the aerodynamic parameters are best obtained from 
experimental section data for the lift and pitching moment of the given 
section at the Mach number of interest. If such data are not available, 
the section lift— curve slope m Q may be approximated by the G-lauert— 
Prandtl and Ackeret relations in the subsonic and supersonic regions, 
respectively, and by an arbitrary constant value which depends on the 
critical Mach number in the transonic region, as indicated' in figure 6. 

The aerodynamic center is located near the quarter— chord point at subsonic 
speeds and at 4-0 to 4-5 percent of the chord at supersonic speeds. Its 
location at transonic speeds Is not well established; it may move forward 
of the quarter— chord point, however, before receding to the supersonic 
position. 



8 


NACA TN No. 1680 


Use of the Curves 

The aerodynamic parameters are selected, for a given Mach n umb er and. 
the section lift— curve slope is corrected for finite- span effects in the 
manner previously indicated. With these parameters and the given geometric 
parameters the ratio dp/ap may he determined from equation (3)* The 

value of ap or dp is then obtained from figure 4 and the divergence 
dynamic pressure calculated from either of the two following relations: 


0.D = a D 


(G-J)r L 

cos 2 A e l c r 


= d T 


(Bl)r 

mgCyL^ cos 2 A 


cot A 


The divergence speed is given by the relation 


= 


i2_ 

p/2 


(5) 

(5a) 


(6) 


If the value, of qp as calculated in this manner is negative the 
wing cannot diverge, since a negative dynamic pressure does not correspond 
to any real speed. A negative value, however, is still useful as a 
reference value. 

Unless the wing under consideration is solid or has geometrically 
similar cross sections along the span, it must be kept in mind that the 
actual divergence speed or dynamic pressure may be below that calculated 
by this method for reasons cited subsequently in the "DISCUSSIOU.” 


Use of Approximate Formulas 


The relation between the values ap and dp may be approximated 
closely by a straight line. The agreement between the exact curves and 
curves representing the linear relation is seen in figures 2 and 3 in the 
case of a uniform wing. From the linear approximation an expression for 
the divergence dynamic pressure may be obtained of the following form: 


<lD = 


(GJ); 


m e c r L3 


cos 






1 - S 2 


(Gj) r / L ' 

(EI) r \®i c r> 


tan A 


(T) 


T-T 
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where the constants Kq and Xq are given in the following table: 


X 

% 

*2 

0.2 

2.8l 

o.6i4 

..5 

2.74 

.497 

1.0 

2.47 

.390 

1.5 

2.22 

.326 


Equation (7) may he used with almost the same accuracy as the curves 
of figure This equation should he particularly convenient for calculating 
the stiffnesses required for a given value of the dynamic pressure at 
divergence. 


EXH52UMEITIAL RESULTS 


As a check on the theory several divergence teBts on uniform swept 
wings were made in the Langley 4. 5— foot flutter research tunnel. In the 
first series of tests a thin plate of 24S— T aluminum alloy, 5 inches 
bv 30 inches hy 0.126 inch, was held in the tunnel essentially as 
shown In figure l(a); the angle of sweep was varied hy means of the 
rotating root fixture. The experimental divergence dynamic pressure was 
taken to he the highest value at which the wing would remain In an 
undeflected position when the root section was at zero angle of attack 
relative to the true air stream. From the experimental data obtained, 
shown in table I, the divergence dynamic pressure is plotted against the 
angle of sweep in figure 7(a). The variation of divergence dynamic 
pressure with sweep obtained from the theory of the present paper Is 
also shown in figure 7 (a) . No experimental data have been obtained for 
sweephack; the wing fluttered, rather than diverged, at zero sweep, and 
flutter would undoubtedly be critical at any angle of sweepback. However, 
the agreement between the present theory and experiment for angles of 
sweepforward up to 4-0° is excellent. 

A similar series of tests was run for another rectangular al u minum 
wing, 4 Inches by 24 inches by 0.0977 inch. (See table I.) The results 
given In figure 7(b) also show good agreement between theory and experiment. 

A third series of tests was run on flat— plate aluminum models, 

5 Inches by 30 inches by 0.125 inch, that were not rectangular but were 
clamped at the root and cut at the tip parallel to the wind stream 
(fig. 1(c)). A separate model was used for each angle of sweep, but 
the length along the leading edge was held constant at 30 inches . These 
wing plan forms do not correspond to those assumed for the theory, but 
they do represent more closely conventional swept-wing plan forms. The 
experimental results are presented in table I and plotted in figure 7(c ) . 
The theoretical variation of the divergence dynamic pressure with sweep 
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obtained by using the effective root discussed previously ( rc -- . . 
is shown, by the solid— line curve. The variation obtained frcn • ■ . 
the effective root is considered to be at the line through the Us.- -r 
point of the leading edge (root B of fig. 7(c)) is shown by the- 
curve’. It appears that the use of root A yields conservative rejul-.j 
all sweep angles of practical concern, but that the’ results obtain- i L - 
using root B are on the average in somewhat better agreement viuh -,h. ; 
data for low ^ind medium, angles of sveepforward. 

The curves of figure 7 are based on calculated values of (Gj)v. 
and (El) r . The section aerodynamic center was assumed to be at 
the c j\ point, so that ep is 0.25. A section lift slope of 2ic vau 

assumed and corrected for span effects as shown in equation (ha) . 
Compressibility effects were neglected in all comparisons, since the 
test Mach numbers were too low to warrant a correction. 


■ } 


i 


P 






DISCOSSIOU 


The aerodynamic and structural assumptions made impose limitations 
on the accuracy of the analysis. The aerodynamic assumptions are 
concerned with the magnitude of the deflections involved as well as the 
treatment of the aerodynamic induction and compressibility effects. The 
assumption of small deflections made in obtainipg the effective velocity 
component (see appendix, equation (A5a)) yields results which are 

too low at high angles of sweep. The air forces are consequently 
underestimated and the divergence speed is overestimated. The correction 
for aerodynamic induction is only approximate; it is simple to make, 
however, and yields good experimental agreement. For most airplanes the 
divergence speed lies in the transonic or supersonic region, where the 
induction effects are greatly diminished, so that no correction need 
ordinarily be applied. The manner in which compressibility is taken into 
account in the analysis has not been checked experimentally and is 
consequently somewhat uncertain. 



/ 


The location of the aerodynamic center may not he known very 
accurately; fortunately, however, the divergence speed of a swept wing 
is not very sensitive to the aerodynamic— center location at large angles 
of sweep, as may be seen from figure 4(c). The change in dp, and hence 
the change in divergence speed, will be found to be fairly small even for 
rather large changes in the value of ep produced hy changes in the 


aerodynamic-center location or elastic— axis location. At low angles of 
sweep the effect of the location of the aerodynamic center relative to 
the elastic axis is much more pronounced; movement of the elastic axis 
forward or the aerodynamic center rearward will tend to 'raise the 
divergence speed. If the aerodynamic center is behind the elastic axis, 
negative values will usually be obtained for the dynamic pressure at 
divergence, so that divergence will be impossible. 
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The most Important of the structural assumptions involves the existence 
of a straight elastic axis* which permits tending and torsion to be treated 
independently . This assumption appears warranted inmost cases. Although 
the exact location of this axis is often difficult to ascertain, only an 
approximate value is needed in view of the foregoing consideration,, at 
least for large angles of sweep. Hear the wing root, where the elastic- 
axis location is most uncertain, swept wings are so stiff that the location 
is immaterial. 

f 

The assumption of the effective root at the location indicated in 
figure 1(c) (root A of fig. 7( c )) leads to conservative results. This 
assumption appears, therefore, to be preferable to the assumption of the 
root at a point farther outboard, such as root B of figure 7(c). On the 
other hand, the assumption of the effective root B leads to results which, 
at least in the case of the flat— plate models tested, agree somewhat more 
closely with experiments for angles of sweepforward up to about 45°. 

The insensitivity of the divergence speed to the relative location 
of the aerodynamic center and the elastic axis indicates that the effects 
of bending predominate at large angles of sweep (actually at large values 
of dp/ap) . Even at low values of sweep, bending is quite significant; 
the divergence phenomenon of swept wings should therefore be referred 
to as "bending-torsion" divergence rather than "torsional" divergence 
as in the case of unswept wings. 

The stiffness variation used in the analysis for tapered wings 
(El and GJ vary as c^) is realized for wings with geometrically 
similar cross sections; it is obtained for solid wings and closely 
approximates that of actual wings with a taper ratio of the order of Q.2. 

For higher taper ratios, actual wings are more flexible at the tip than 
the fourth-power variation would dictate. In order to investigate this 
effect, the parameter ap has been plotted against taper ratio in figure 8 

for unswept wings with stiffness variations dictated by constant bending- 
stress levels. The assumption has been made that the torsional stiffness 
is proportional to the flexural stiffness which, in turn, is based on a 
load distribution given by strip theory. The computations for the 
divergence speed were performed by a method similar to that of reference 2. 
Figure 8 indicates that the divergence speed of wings with this type of. 
stiffness distribution is in general lower than that obtained for wings 
the stiffness of which follows the quart ic variation. The divergence speed 
of actual wings may be expected to lie between the two curves of figure 8 
and will in general be lower than that for wings the stiffness of which 
follows the fourth-power variation. This point must be considered when 
figure 4 is used; the results furnished by the curves of figure 4 may be 
somewhat unconservative for’ actual wings . 

The effect of sweep may be seen from figures 3 , 4, and 7. Figures 3 
and 4 show that as the sweep angle (and hence the parameter dp/ap) is 
increased negatively, that is, toward increased sweepforward, the 
parameter ap (and hence the divergence speed) decreases rapidly, all 
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other physical parameters remai n ing the same. Beyond a certain angle 'of 
sveepforvard the cos A term in the parameter ap "becomes dominant 

and tends to increase the divergence speed again as seen in figure 7* 
Sveptback vings vith a value of dp/ap larger than about 2 and positive 
values of ep, vhich normally correspond to subsonic speeds, cannot 
diverge; nor can any sveptback or unsvept ving diverge for negative values 
of ep vhich may exist at supersonic speeds. In either case a negative 
divergence dynamic pressure is obtained. This fact does not preclude 
the possibility of flutter in these cases, hovever. 

The effect of taper on the -divergence of svept vings of the assumed 
stiffness variation (GJ and El varying as cM is seen from figure 4. 

For positive values of ep (subsonic case) and positive or small negative 

values of the sveep angle, conventional taper increases the divergence 
speed; for all other configurations it decreases the divergence speed. 

The effects of inverse taper are the opposite of those of conventional 
taper. These considerations must be modified in the caBe of actual 
vings, because of their deviation from the assumed stiffness distribution, 
so that either the curves of figure 4 must be used vith some degree of 
conservatism or more refined analyses must be resorted to. 

Although the lav-speed vind— tunnel tests for divergence have been 
performed only on uniform vings, they serve to corroborate the theoretical 
analysis for uniform vings and hence, indirectly, that for tapered vings, 
since the assumptions are the same “in either case. The rapid decrease 
of divergence speed vith increase in sveepforvard agrees vith the predicted 
variation both qualitatively and quantitatively (up to about 40° sveep- 1 
forvard) . This agreement indicates that the assumptions made concerning 
the structural and aerodynamic behavior (at lav speeds) are justified. 

At values of sveepforvard above about 40° the observed increase in 
divergence speed falls short of the predicted increase. This observation 
agrees qualitatively vith the statement made previously that the analysis 
vould tend to overestimate the divergence speed for large angles of sveep. 

It is difficult to make a direct comparison vith the results of the 
simplified approach of reference 5, but it is quite apparent that since 
the method of this report takes account of the spanvise variation of 
bending and tvisting It should furnish more reliable results. If the 
first natural bending and torsiop. frequencies of a uniform cantilever 

tXD„ 

beam are substituted in the exoression for — = of reference 5. a 

relation may be obtained for qp vhich vill have both the same form as 
equation ( 7 ) and the same value of if a proper finite— span correction 

is applied; the value of Ibj, hovever, vill be 0.299 instead of 0 . 390 ^ 

as determined by the theory of this paper for the unif orm ving (X = 1.0), 

30 that the effect of sveep vill be underestimated somevhat. This 
discrepancy indicates that. If the results of a semirigid analysis are 
applied to actual svept vings, a certain amount of caution must be exercised. 
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CONCLUDIHGr HEMAKKS 


On the "basis of certain assumptions theoretical results hare been 
obtained for the divergence of uniform, and tapered -wings . These results 
are presented in the form of charts and approximate formulas suitable 
for obtaining quick estimates of the divergence dynamic pressure or speed. 
A limited number of low— speed wind— tunnel tests on untapered models give 
good agreement with the theoretical results for angles of sveepforward 
up to 40° And appear to justify the assumptions made in the analysis. 

The results indicate the following conclusions: 

1. The divergence speed drops rapidly as sweepforward increases up 
to about 40°. Wings with sweepback beyond a fairly low value cannot 
diverge . 

2. Moving the elastic axis forward raises the divergence speed 
appreciably for low angles of sweep but' has less effect at higher sweep 
angles. Wings with the elastic axis forward of the aerodynamic center 
(supersonic case) can only diverge for moderate or large angles of 
sweepforward. 

3. For most practical cases, the effect of conventional taper is to 
increase the divergence speed of essentially unswept wings and to decrease 
the divergence speed of wings with moderate and large angles of sweep if 
the stiffness varies as the chord to the fourth power. For stiffness 
variations closer to the ones which are obtained for actual wings, this 
effect may not be observed. In order to obtain accurate results in these 
cases a more refined analysis must be resorted to. 


Langley Memorial Aeronautical Laboratory 

Rational Advisory Committee for Aeronautics 
Langley Field, Va., April l6, 1948 
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APPENDIX 

ANALYSIS OF THE DIVERGENCE OE SWEPT WINGS 


f ' Assumptions 

The analysis inv olves the following limitations and assumptions: 

(a) Aerody nami c induction is taken into account only insofar as 
an over— all correction is applied to strip theory 

(h) Aerodynamic as well as elastic 'forces are based on the 
assumption of small deflections 

(c) The wing is clamped at the root perpendicular to a straight 
elastic axis (see fig. l), and all deformations are 
considered to he given by the elementary theories of 
bending and of torsion about the elastic axis 


Aerodynamic Forces 

In keeping with assun^t ions (a) and (b) , the force per unit width 
on a wing section in a plane perpendicular to the elastic axiB resulting 
from bending and twisting deformation! is given by 

l =m e a Q PT A 2 c (Al) 

. where ou and V* are the effective angle of attack and the component 

^ IV 

of the free— stream velocity in the plane perpendicular to the elastic 
axis, respectively. If the free— stream velocity is resolved into three 
components, one along the local tangent to the elastic axis , one parallel 
to the chord, and one perpendicular to the other two, the effective 
angle of attack may be obtained as the ratio of the third component to 
the second and the effective velocity as the vector sum of the two 
components . Thus , 

sin cp cos A — cos 9 sin T sin A . 

tan a e = (A 2 ) 

cos cp cos A + sin 9 sin T sin A 

and 

Yj^ = Y 2 j^(cos 9 cos A + sin 9 sin T sin A) 2 

+ (sin 9 cos A — cos 9 sin T sin A)^j (A 3 ) 

The most convenient way of obtaining these relations is probably by use 
of vector analysis. 
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For small values of 9 and T these relations simplify to 

a e =9 - r tan A (A 4 ) 

T a = V cos A /l + a e 2 (A 5 ) 

= V cos A (A 5 a) 

so thAt, for small deflections, equation (Al) can. he expressed as follows: 

2 = (cp — T tan A) m e £'V 2 co8%c (a6) 


The moment of this force about the elastic axis is assumed to he 
given hy the product Ze^c. The presence of a moment at zero lift due 

to camber will not affect the results in any manner, since the divergence 
speed depends only on the rates at which the aerodynamic and structural 
forces increase with the deformations. 


Differential Equations 

The differential equations of bending and torsion of the wing 
referred to the elastic axis are 


dy 2 

_d_ 

dy 


( EI §) - 1 
( v v) ‘ - 


(AT) 


(AS) 


respectively. Substituting the smai 1 -deflection expression for 2 
(equation (A6)) into equations (AT) and (A8) yields the two simultaneous 
differential equations: 


_df_ 

dy 2 


(A9) 


^EI = ^Y^os^AEiigC (cp — T tan A) 


— (crj =— ^V^cos 2 Am e-,c^ (9 — T tan A) (AlO) 

dy \ dy/ 2 e x • 


These equations are subject to the following boundary conditions: 


Zero local dihedral and twist at the root, 

r (0) =0 
9 (0) = 0 


(All) 

(Alla) 
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Zero moment, torque, and. shear at the free end. 


El ^ 
V *3/ 


= 0 


(A12) 


y=L- 


GJ ^2. 


= 0 


(A12a) 


y=D 


EI 


a 2 r 

dy ? 


= 0 


(A12b) 


y=L 


Equations (A9) and (A10) may he solved hy numerical methods for 
any arhitrary stiffness and chord variation. In the case of untapered 
wings and linearly tapered -wings with the stiffnesses varying as the 
chord to the fourth power, they can he solved directly, as shown in the ■ 
following sections. 


Solution for Uniform Wings 

* 

If the tending stiffness, the torsional stiffness, and the chord 
of the wing have constant values of (El) r , (GJ) r , and c r , respectively, 
along the wing span, equations (A9) and (A10) become 


r ,,t tan A = d(cp — T tan A) 

■ cp *’ = — a(cp — F tan A) 


(A13) 

(AUO 


where the differentiation denoted hy the prime is with respect to 
and the two dimensionless parameters a and d are defined hy 

q cos^AniQeqc^L 2 

a= m; 


T) = 


(A15) 


q cos^mgCpL^an A 

d = 

(EI) r 


(A16) 


Differentiating equation (Al4) once and combining it with equation (A13) 
yields the single differential equation 

Og’ 1 1 + aa e ! + &a e = 0 (A17) 


where — r tan A. 


t-»M 
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The general solution of this equation is 


% = /_ A i© 

i=i 

where the r^*s ‘are the roots of the characteristic equation 

r^ + ar + d = 0 

and the A^’s are arbitrary constants. 

The boundary conditions for equation (A17) are: 

From equations (All) and (Alla) 

Oe(0) = 0 

From, equations (A12) and (A12a) 

V (1) = 0 

From equations (AlA) and (A12b) 

1 (l) + aog(l) =0 

Subs + ituting equation (Al8) into these boundary conditions gives 

^A i = 0 

i=l 

' r 5 - r,.A,e rl -0 


(A18) 


' (A19) 


(A20) 


(A20a) 


(A20b) 


’ £~ T,A,e ri . 0 
i=l 

(rj_ 2 + a.)k±e T± = 0 


(A?l) 


The condition for divergence is that there be a nonvanishing 
solution for the A^sj that is, a solution for which <Xq is not 

zero along the entire Bpan. Therefore, one or more of the Aj_*s must 
be different from zero (see equation (Al8)). Hence, the deter mina nt of 
the coefficients of the Aj_ l s in equations (A21) must vanish. Thus, 
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1 


l 


1 


r l e 


r l 


r 2 e 


r 2 


r^e r 3 


= 0 (A22) 


+ e 


a )' 


r l 


(r 2 2 + a)e r2 (r-^ 2 + a)e r 3 


Critipal combinations ap and d D of the parameters a and d 
are obtained if the combinations give rise to roots r^ r 2 , and r^ 

for equation (A19) such that equation (A22) or its expanded equivalent 
is satisfied. The curve of ap . against dp in figure 2 constitutes a 
plot of these critical comb ina tions ; points on the curve were computed 
by assuming values of ap and solving for djj by trial. 


Solution for Tapered Wings 

For tapered -wings the chord varies linearly and the bending and 
torsion stiffnesses are assumed to vary as the fourth power of the chord. 
Thenj 

c = kc r 


where 

k = 1 - (1 - X)i] 

and X is the taper ratio c-fcAr • Furthermore } 

GJ = (Gj) r k^ 

EE = (EI)^ 


The differential equations (A9) and (A10) then become 
tan A[k 3 r m + Sh 2 !’ 11 + 12kr* - dr^r] + dqcp = 0 
|k 2 <p n + 43sp* + aijtpJ - aijiT tan A = 0 

where 


Sjji — 


a 

(1 - x) 2 


dp — 


d 

(1 - X)3 


(A23) 

{k2h) 
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The differentiations denoted by the primes in equations (A23) and (A24) 
are with respect to k rather than t}. This .procedure places the 
differential equations in the form of the Euler (or Cauchy) equations , 
which are easily tractable mathematically . 


Differentiating equation (A24) once., multiplying it by k, and 
combining the result with equations (A23) and (A2.k) yields a single 
differential equation in <x e : 

-k^ag" 1 + 8k 2 Qg lJ + (12 + sqOkctQ* + (28^ — dqOctg = 0 (A25) 


The solution of this equation is 


CL = 

e 



(A26) 


where the s^’s are the roots of 


s(s — l)(s — 2) + 8s(s — l) + (12 + aqr)s + (2aj — dtj) = 0 


or 


S^ + 5s 2 + (6 + 8rji) s + (2aqi — dqr) = 0 

and the Bi*s are arbitrary constants.. 

The boundary conditions are: 

From equations (All) and (Alla) 

a e (l) = 0 

From equations (All) and (A12a) 

ct e T (X) = 0 

From equations (A24), (A12a) and (A12b) 

(X) + s^gCX) = 0 

Substituting equation (A26) into the boundary conditions gives 


i=l 

3 


s^X® 1 = 0 


i=l 


IlL [ s ± C s i - 1) + Sqil Bj_X Sl = 0 

-? *- -L -l 


i=l 


(A2T) 


(A28) 

(A28a) 

(A28b) 


(A29) 
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Setting the determinant of the coefficients of tlie B.*s in 
equation (A29) equal to zero yields, after some simplification. 


1 

1 

1 



x s i 


s 3 X S 3 

= 0 

(A30) 

+ a^X 8 ! 

(s^ + aji^X S 2 

(s^ 2 + ap)x S 3 




Equation (A30) is very similar to the corresponding determinants! 
equation (A22) for untapered wings. For particular values of the taper 
ratio 1, it determines -critical combinations of a and d. Calcu- 
lations have been carried out for taper ratios of 0.2, 0.5, and 1.5j sna- 
the results are given by the curves of ap ' against dp/ap and dp 
against ap/dp in figures 4(a), Mb), and 4(c), together with the 
results of the analysis for the uniform wing (taper ratio 1.0). 
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TABEE I.- RESULTS OF DIYERGERCE TESTS 


-A 

(deg) 


(experiment) 
(lh/sq ft) 


(calculated) 
(lh/sq ft) 


Series 1 


5 in. "by 30 in. By 0.126 in.; square tips; 0! = 0.25; 
(GJ) r = 13,330 lB-in. 2 ; (El) r = 8,830 It^in. 2 



a 112.80 

83.20 

44.21 
26.50 
24. 4l 
24.84 
24.95 
26.65 


178.6 

80.7 
40.9 

27.0 
26.2 

31.1 
31.1 

39.7 


Series 2 


4 in. By 24 in. By 0.0977 in.; square tips; = 0.25; 
(GJ) r = 4,980 IB— in. 2 ; (El) r = 3,300 lB-in . 2 


0 

a l06.90 

163.1 

5.0 

76.99 

73.8 

14.7 

43.10 

37.4 

30.0 

26.51 

24.7 

45.0 

22.44 

23.8 

60.0 

23.36 

32.0 

69.6 

' 23.74 

50.5 

Series 3 


5 in. By 30 in. By 0.125 in.; tips parallel to plane 
of symmetry; e-]_ = 0.25; 

(Gj) r = 13,050 IB— in. 2 ; (El) r = 8,560 lB-in. 2 ' 


5.6 

15.5 
30.1 ‘ 
45.2 
3 59.7 

74.6 



79 

4l 

31 

28 

38 


54.62 


Root A 

Root B 

74.0 

38.3 
26.5 1 

25.3 1 
33.7 
79-3 

74.8 

40.8 
30.3 

32.8 
53.1 

232.4 


^Flutter . 

^Model thickness, 0.123 inch. 

G Corrected for discreoancy in thickness 
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Figure 2.- The relation between the parameters and 
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